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Abstract. A change of the component separation in the profiles of the 
binary pulsar PSR B 1913+ 16 has been observed for the first time (Kramer 
1998) as expected by geodetic precession. In this work we extend the 
previous work by accounting for recent data from the Effelsberg 100- 
m telescope and Arecibo Observatory and testing model predictions. We 
demonstrate how the new information will provide additional information 
on the solutions of the system geometry. 



1. Introduction 

The binary pulsar PSR B1913+16 exhibits a measurable amount of geodetic 
precession due to a misalignment of its spin axis and orbital angular momentum 
(relativistic spin orbit coupling) by an angle 5 (see also Stairs et al., for geodetic 
precession in PSR B1534+12). The angle between the line of sight and the spin 
axis, C, changes according to: 

cos C(t) = — cos 6 cos i + sin 5 sin i cos f2 prec (t — To) , 

where i is the orbital inclination of the binary system. General Relativity (GR) 
predicts a rate of precession f2 pre c = 1-21 deg yr _1 , leading to a change in the 
pulse profile and polarization characteristics with time. 

Previous observations by Weisberg, Romani and Taylor (1989) detected a change 
in the relative amplitude of the trailing and leading component, but no change 
in the component separation. Also, there were no changes in the position angle 
swing of the linear polarization noted (Cordes, Wasserman & Blaskiewicz 1990). 



2. Recent Observations Cast Some Light 

The evolution of the observed radio pulse profile of PSR B1913+16 would be 
determined when the four parameters a, 5, p, Tq were specified for a given i, where 
a is the angle between the magnetic pole and the spin axis, 5 is as mentioned 
above, p is the angular radius of the emission cone and To the epoch of precession. 
Based on recent data from the Effelsberg 100-m telescope, Kramer (1998) first 
noticed a change in the component separation W. This allows a solution region 
in parametric space that respects GR and the hollow cone model (HCM). 



W(t) = 2arccos 



cos p — cos a cos ((t) 



sinasin£(i) 
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Effelsberg polarization data has been used in respect with the rotating vector 
model (RVM), to further narrow down the solution region. 



ip(t) = arctan 



sin a sin 



cos a sin — sin a cos ((t) 



+ const. 



The inclination of the binary orbit with respect to the line-of-sight, i, is 



known from timing observations modulo the ambiguity i 
possibilities are: Case A: i = 47.2°, Case B: i = 132.8° 
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3. The Solution Today 
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Figure 1. The pulse component sep- 
aration data (Effelsberg + Arecibo) in 
conjunction with the HCM was used to 
produce a \ 2 area plot indicating the 
3cr confidence level solutions for a and 
S, areas inside the surfaces traced by the 
black lines. All possible solutions were 
then fed into the RVM to produce posi- 
tion angle curves, which were compared 
to Effelsberg polarization data, only al- 
lowing for solutions in the unshaded ar- 
eas. So, the solution must be located 
on a black surface inside an unshaded 
region (for detailed explanation of 
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The set of plausible solutions has been well limited, as shown in Fig 1. 
The inclination of the binary orbit is most likely to be Case B rather than 
Case A (Wex, Kalogera, Kramer 2000). More precise use of the polarization 
data plus a molding of the two models (HCM & RVM) to specifically suit PSR 
B1913+16 (see also Weisberg & Taylor in these proceedings) could further limit 
the solutions and enable a prediction on the future of this system. 
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